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Abstract: The realization of global symmetries can depend on the geometry of the un- 
derlying space. In particular, compactification can lead to spontaneous breaking of such 
symmetries. Four-dimensional QCD with fundamental representation fermions embedded 
in a space with one compact spatial dimension has a critical length, at which the theory 
undergoes a phase transition and develops a ground state that is no longer charge conju- 
gation invariant. We show this behavior with simulations of three color, four flavor QCD. 
We use unrooted staggered fermion at two values of the lattice spacing and several quark 
masses. We discuss the dependence of the transition on the dynamical fermion mass as 
well as its connection to the finite temperature and chiral phase transitions. 
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1. Introduction 
1.1 Motivation 

When field theories are defined on spaces with one or more compact directions, the sizes 
of the compact dimensions can influence the existence and location of phase transitions. 
The most familiar example of this behavior is the use of a compact temporal dimension to 
study a field theory at finite temperature. In this paper we describe such a situation, which 
occurs in ordinary QCD (with dynamical fermions): a sufficiently small compact spatial 
dimension induces a transition into a phase in which charge conjugation is spontaneously 
broken. The existence of this new phase, and of the length scale which characterizes the 
transition, depend on the choice of periodic fermionic boundary conditions in the compact 
dimension(s). 

The effect we describe is similar to the deconfinement transition in QCD with dy- 
namical fermions. The action of SU (N) pure gauge theories is invariant under gauge 
transformations which are periodic up to multiplication by the center of the gauge group, 
Z{N) for SU(N). These theories have an order parameter, the Polyakov loop or Wilson 
line, which is sensitive to such transformations. The low temperature phase of the theory 
is one in which the vacuum expectation value of the order parameter vanishes. In the 
high temperature phase, the order parameter spontaneously magnetizes along one of the 
directions of the center of the gauge group, signaling deconfinement. Now add dynam- 
ical fermions. The Euclidean path integral for fermions has a probability interpretation 
when the temporal boundary conditions of the fermions are chosen to be anti-periodic. 
The fermion action (for fundamental representation fermions) explicitly breaks the Z(N) 
symmetry, favoring the vacuum in which the Polyakov loop takes a positive value. 
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For the physical case of three colors, the deconfinement transition is first order for the 
pure gauge theory and is believed to disappear as the fermion mass is reduced from infinity, 
as only one vacuum configuration is energetically favored. This is all well known: for an 
early review, see [1]. 

However, if the fermions obey periodic boundary conditions, the vacuum at the positive 
center element is disfavored. If the number of colors is odd, there are two degenerate minima 
in the free energy, corresponding to complex-conjugate values of the order parameter. One 
of these is selected as the vacuum and thus charge conjugation is spontaneously broken. The 
radius of the compact dimension at which the transition from C-conserving to C-broken 
occurs is a physical property of QCD. 

The earliest discussion of this new phase that we are aware of is by van Baal [2,3]. He 
was interested in QCD with all three spatial dimensions taken to be small. In this case 
there are 2 3 degenerate vacua, two per small dimension. While he suggested that numerical 
simulations of QCD with compact dimensions might be interesting, we are aware of no 
follow-up studies. 

Our interest in this problem was initiated from a quite different source: a discussion of 
the validity of the "orientifold QCD" programme of Armoni, Shifman, and Veneziano [4-6]. 
These authors proposed a large- iV c equivalence between 7V=1 super Yang-Mills and QCD 
with a single fermion flavor in the symmetric or antisymmetric tensor representation. The 
latter theory is equivalent to one-flavor QCD for N c = 3. The authors used this connection 
to make a successful prediction of the value of the condensate of one-flavor QCD [6,7]. 

This summer, Unsal and Yaffe [8] pointed out that a necessary condition for the non- 
perturbative equivalence of different theories is that global symmetries are realized in the 
same manner in both. They show that if the theories are considered on IR 3 x S 1 and the 
radius of S 1 is chosen sufficiently small, QCD with tensor representation fermions will 
spontaneously break charge conjugation symmetry, while the supersymmetric theory does 
not. This work is a specific application of their earlier, more general studies [9-11]. 

Their discussion is very general: it depends only on the gauge group and on the rep- 
resentation of the fermions. We decided to look for the C-violating phase in numerical 
simulations, and to do this for a choice of parameters which was computationally the most 
inexpensive. Accordingly, we performed simulations with three colors and four flavors of 
fundamental representation dynamical fermions. We concentrate on zero-temperature sim- 
ulations with one compact dimension (this is done in the simulation by taking the compact 
dimension much smaller than the other ones) but also briefly consider other situations: 
finite temperature plus one or more compact spatial dimensions. The transition is strong 
and easy to observe. We have not attempted to determine its order. 

In the remainder of this section we outline theoretical expectations for C-odd phases, 
and then proceed to simulations and their results. 

1.2 Weak and strong coupling 

Asymptotic freedom ensures that a perturbative calculation of the effective potential for 
Polyakov lines becomes reliable if their length L is short enough. In this case the coupling 
at the scale 1/L is small and higher order corrections become negligible. For nj flavors of 
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massless Dirac fermions in the fundamental representation of SU(N) the effective potential 
for the Polyakov loop P is given by [8] 
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where D 2 is the covariant Laplacian for a constant background field in a given represen- 
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where [v] = t!mod2-7r. In the case of periodic boundary conditions the effective potential 
has a unique minimum at [v] = tt, while thermal boundary conditions favor a Polyakov loop 
on the positive real axis. Gauge invariance then dictates that the actual quantum vacuum 
configurations are those center elements of SU(N) closest to v = and v = ir, respectively. 
We thus expect a spontaneous breaking of charge conjugation for SU(N) with N = 2k + 1, 
k > 1. 

A few comments are in order 1 : For gauge group U(N) the phases v = and v = it 
are elements of the center and thus there is no spontaneous breaking of charge conjugation 
for fundamental representation fermions. This is different from the case of tensor repre- 
sentation fermions, which polarize the Polyakov line to a phase of ±7r and thus always 
break charge conjugation. Moreover, for large N the vacua for fundamental representation 
fermions and periodic boundary conditions are separated from ir only by 1/N effects. Thus 
the symmetry breaking disappears in the large ./V limit together with the difference between 
U(N) and SU(N) gauge groups and the effects of fundamental representation fermions in 
general. 

Another calculation, by Hollowood and Naqvi [12], of the phase structure of gauge 
plus tensor-representation fermion theories with three small, compact spatial dimensions, 
shows the presence of a confinement-deconfinement transition, irrespective of the fermionic 
boundary conditions. 

The variation on this explanation which is probably most familiar to a lattice prac- 
titioner comes from strong coupling and large mass expansions. (A probably incomplete 
set of early references includes [13-20].) Consider a lattice which is long in the thermal 
direction and in all but one of the spatial directions, which we denote by x. There are N x 
lattice spacings in the short direction. We parameterize our fermion action in a way which 
is suitable for a hopping parameter expansion, integrate out the fermions and then expand 
the determinant as a power series in the hopping parameter k. At large mass, k ~ e ~ ma 



*We thank Mithat Unsal for pointing those out to us. 
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where a is the lattice spacing. In lowest order, we find an effective action whose leading 
term is proportional to 



where P (re) is the Polyakov line oriented in direction x and n labels all sites in the lattice 



the fermions in direction x: negative for anti-periodic, positive for periodic. This result is 
qualitatively similar to the result of the calculation of the effective potential: the effect of 
a short direction with anti-periodic boundary conditions is to break the center symmetry 
to its trivial value, while periodic boundary conditions bias the Polyakov loop toward a 
negative value. If —1 is an element of the center, then there is again a unique minimum; 
otherwise, there will be degenerate vacua. 

We conclude this section with a few simple and general remarks: 

• The actual location of a phase transition (the value of the critical compact radius L c ) 
can only be determined by a non-perturbative analysis. Its value is a fundamental 
property of QCD, like a hadron mass. 

• As the quark mass goes to infinity, L c = 1/T C , the critical temperature for deconfine- 
ment in the quenched theory. 

• As m q — > there is presumably a chiral symmetry restoration transition at some 
L c . All of the old questions of the interplay of deconfinement and chiral symmetry 
breaking, such as whether there are two separate transitions, or one, are present 
again. Unsal and Yaffe have pictures of proposed phase diagrams with separated 
transitions. 

• Finally, for SU(3), the pure gauge deconfinement transition is first order. Fermions 
break the Z(3) symmetry to Z(2). A first order transition is stable under small 
symmetry breaking effects, so we expect the transition to persist as the quark mass 
is lowered from infinity. The transition might end in a second order critical point, or 
it might convert to a line of second order points. In either case, we expect that the 
second-order transition would be in the universality class of the three-dimensional 
Ising model. 

2. Simulations 
2.1 Methodology 

In our simulations we employ the standard setup for dynamical fermion lattice simulations; 
we have three spatial dimensions with periodic boundary conditions and one time dimen- 
sion with anti-periodic boundary conditions for the fermions. The geometry of 1R 3 x S 1 
is approximated by keeping all directions but one of the spatial ones (the x direction) 
large. We check this by verifying that Polyakov lines wrapping the large directions are 
consistent with zero and show no phase preference. Note that the only difference to a finite 
temperature simulation is the choice of boundary conditions in the short direction. 
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with x = 0. The sign of the right hand side of Eq. (1.4) reflects the boundary conditions of 
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Table 1: Simulation parameters. 



In the present study we choose improved (unrooted, i.e. four-taste) staggered fermions 
as a computationally inexpensive way to do simulations with large dynamical fermion 
effects. We work with a state-of-the-art improved action to minimize cutoff effects. (We 
have also performed simulations with unimproved staggered fermions. They show similar 
results for a transition, but are disfavored for estimating continuum numbers.) While in 
principle there are many parameters one could vary (temperature, length and number of 
compact directions, quark mass, flavor number) we will focus on four-flavor simulations 
at zero temperature with one compact spatial dimension and only vary its length and 



the quark mass. Other choices are briefly discussed in Sect. 3JJ . We employ the Hybrid 
Monte Carlo algorithm from the publicly available MILC code 2 for improved staggered 
quarks [21-23] on a Symanzik gauge background. 

The parameters of our simulations are summarized in Table [|, where the value of f3 
refers to the Symanzik improved gauge action. We ran 500-2000 trajectories per coupling 
value for the N x = 4 simulations and 200-500 at N x = 6. In all simulations the tadpole 
coefficient uo is fixed to a value of 0.875 for all simulations instead of determining it 
self-consistently at each value of (3 and am. In the scale setting runs we measure HYP 
smeared [24] Wilson loops to extract the hadronic scale ro [25] from the static quark 
potential [26] with results shown in Fig. |l|. The solid lines are parameterizations of the 
form 

log(r /o) = c + ct/3 + c 2 p 2 , (2.1) 

which we use to interpolate ro- Throughout this work we will assume ro = 0.5 fm. Note 
that since we do not tune uq self-consistently, our bare action differs from the one used in 
Ref. [27] and we therefore obtain rather different values of r$/a as a function of /?. 

2.2 Results 

We first want to establish that we can indeed find a set of parameters where the results from 
Sect. |1.2j hold qualitatively. To this end we make a comparison: we thermalize two lattices. 
One has a compact spatial direction, and our initial gauge configuration was fixed to the 
identity. The second configuration has a compact temporal direction, and its initial gauge 



2 http : //www. physics .Utah. edu/~detar/milc/ 
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Figure 1: The hadronic scale ro/a as a function of (3 in the ranges relevant for the simulations 
listed in Table 0. 




Figure 2: Thermalization histories of the Polyakov line in the short direction. Left panel: short 
time direction. Right panel: short spatial (x) direction. 



fields had a time-like Polyakov loop with a phase of exp(— 2m/3). We choose N t / X =A for 
the compact direction, which corresponds to ~ 0.6 fm at /3 = 6.65 and a quark mass of 
am = 0.2. In both cases the non-compact directions all have 10 lattice sites. 

In Fig. ^ the direction of thermalization is indicated by an arrow for both cases. One 
can clearly see that the finite temperature simulation (left panel) finds a stable minimum 
on the real axis. The accumulation on the complex Z{2>) leaf (indicated by a circle) might 
indicate a meta-stable state. The simulation with compact spatial dimension (right panel) 
shows the opposite behavior. After following the real axis to the origin it settles to an 
equilibrium (indicated by a dot) with a phase around —2ir/3. Thus we encounter precisely 
the behavior that the analytical arguments from Sect. |1.2j suggest and now proceed to a 
more detailed study of the transition. 

We focus on simulations with one compact spatial dimension (x) at zero temperature 
and vary the extent of the compact dimension by changing (3. Fig. ||] shows scatter plots 
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of P x after equilibration from sets C2. For better readability each point represents the 
average of 5 consecutive measurements of P x . The length of the compact dimension is 
indicated in units of tq. The numerical data confirms the expectation that the fermion 
determinant pushes the Polyakov lines toward the negative real axis. At the largest extent 
we simulated with N x = 4, corresponding to L x = 2.48 ro, P x is centered at -0.0275(17) on 
the real axis. With decreasing length, it moves to more negative values and spreads in the 
imaginary direction toward the complex center elements. 

Below a certain critical length, the distribution becomes disjoint and one of the complex 
Z(3) elements is selected as the ground state configuration with only the magnitude of P x 
increasing with decreasing L x . Thus, for SU(3) the fermion polarization spontaneously 
breaks charge conjugation invariance, indicated by a non-vanishing imaginary part of the 
Polyakov loop. With the chosen algorithm and at the volumes we simulated, no tunneling 
events between the two vacua were observed. This situation changes in finite temperature 



simulations and/or with more than one compact dimension, see Sect. 3.1. 

In principle the imaginary part of the Polyakov loop could be used to make quantitative 
statements about the transition. However, close to the critical length the increased spread 
of the Polyakov lines (see e.g. the L x = 1.74 data in Fig. |||) would only affect the variance 
of the imaginary part. An observable better suited to make quantitative statements is the 
correlation of the phase of the Polyakov loop with that of one of the group centers e 2km ^ 3 . 
Another advantage is that, unlike the imaginary part, the phase of the Polyakov line does 
not require renormalization. We map the phase range between two SU(3) center elements 
to the full circle by taking (P/|P|) 3 and then project onto the real axis, 

S{P) = Re(P/\P\) 3 = cos(3argP) . (2.2) 

The expectation value of S(P) is close to one in the broken phase, when all Polyakov 
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lines are in the vicinity of one of the Z{3) centers. It is zero for randomly distributee- 
phases and negative if the phases lie between the Z{3) elements, e.g. around the negative 
real axis. (This observable is reminiscent of a similar one used to locate the pure gauge 
transition [28].) The center correlation S(P X ) for the N x = 4 simulations is shown in Fig. [| 
as a function of the length of the compact dimension. For the dynamical simulations we 
observe that the value of S abruptly decreases from almost unity to small negative values 
before approaching zero in the large volume limit. For the heavier data sets C3 and C2, 
the transition occurs at ~ 0.8 fm and for the lightest set CI with am = 0.05 the critical 
length is slightly larger at 0.93 fm. The extraction of the critical length is discussed in the 
next section. 

For comparison we also show the quenched results, run Q. As discussed above, the 
quenched limit of our setup is just the familiar finite temperature phase transition for the 
pure gauge theory. Here the behavior is qualitatively different as all three center elements 
are equivalent vacua in the broken phase. Moreover, the heatbath/overrelaxation algorithm 
allows tunneling between those. With our chosen observable the transition appears much 
smoother than in the presence of dynamical fermions. 

To establish that we are indeed seeing a physical effect and not just lattice artifacts, 
we simulated finer lattices (N X = G). Data from the N X = Q simulations are shown in Figs. || 
and ||. 

The results from the finer lattice spacing are qualitatively similar to the N x = 4 data 
presented in Figs. ^ and The transition is still present and only slightly more rounded. 
The critical length has moved to somewhat larger values and there is no longer a region of 
negative values of S(P X ) just above the transition. One might thus interpret this behavior 
as a lattice artifact that disappears when going from N x =4 to 6. 

Finally, we consider the chiral condensate (ijjip)- A non-vanishing expectation value 
indicates spontaneous breaking of chiral symmetry, which occurs in large volumes. Finite 
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Table 2: Results for the critical length. 



volume, non-zero quark mass as well as discretization errors tend to wash out this transition 
such that only a smooth crossover is observed. As we are only interested in the qualitative 
behavior across the transition, we consider the unrenormalized lattice condensate. The 
results from the N X =A runs, again as a function of L/ro, are shown in Fig. [7|. As expected, 
for the heavier masses (C2, C3) no strong signal is observed. However, as marked by 
the lines, there is an inflection point, where the slope (as a function of L x ) increases as 
the compact dimension gets bigger. The former is located just below L = 1.5ro, which is 
where the center correlation (see Fig. |j) starts to drop. The lightest mass (Cl) shows a 
much stronger crossover and again the inflection point agrees with the data from the center 
correlation. Almost identical results are obtained at the finer lattices. 

Thus the data indicate that in Nf = 4 QCD, at 
zero temperature with one compact spatial dimen- 
sion, there is a single crossover, i.e. a critical length, 
below which chiral symmetry is restored and charge 
conjugation is spontaneously broken. 
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3. Interpretation 



Figure 7: The scalar condensate 

vs. the extent of the compact dimension 

for the N r = 4 runs. 



As a definition of the critical coupling j3 c we take 
the value, where the center correlation S'(P) = 0.5. 
In practice we linearly interpolate between the two 
values of j3 above and below this point. Since this 
quantity shows a sharp transition the results will 
not change significantly when different definitions 
are adopted. The values of (3 C and the critical length 
L c for all runs are given in Table @. The error on L c includes the (statistically independent) 
error of ro- 

Taking into account the rather coarse lattice spacing and the arbitrariness in the precise 
definition of the critical length, the result for the quenched case is in good agreement 
with other determinations of the quenched critical temperature (see e.g. [29], who quotes 
T c ~ 270 MeV.) Comparing the data from the coarse dynamical lattices in Table || we see 
that the critical length increases with decreasing quark mass. Data at fixed N x -am = mL c , 
i.e. Cl/Fl and C2/F2, suggest that scaling violations are not large and that the lattice 
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0.6 0.7 0.8 0.9 

L c /fm 

Figure 8: The phase diagram at N x — 
4, where we have three dynamical and a 
quenched data point for the critical length 
of the compact dimension. 




Figure 9: Equilibrated Polyakov lines for 
13 = 6.65, am = 0.2 and [N x ,N y ,N z ,N t ] = 
[4,10,10,4]. 



artifacts decrease the value of L c . An extrapolation assuming a 2 scaling gives continuum 
values of L c = 1.027(26) fm and 0.922(29) fm for mL c = 0.2 and 0.8, respectively. 

Using all the N x = 4 data, we obtain a phase diagram of the critical length as a 
function of the quark mass as shown in Fig. ||. The critical length increases as the quark 
mass is lowered from infinity and although we have no data at very light quark masses, the 
dependence does not seem to be very strong. Clearly, more quantitative statements will 
require simulations at finer lattice spacing and more (and lighter) quark masses. 

3.1 Other geometries 

To briefly explore other geometries, we chose a parameter set in the broken phase {am = 0.2, 
(3 = 6.65) and simulated all possible combinations of compact (N = 4) and non-compact 
(N = 10) dimensions, including finite temperature. In all cases we confirm our expectation 
that compact spatial dimensions support a Polyakov loop in one of the two complex Z(3) 
directions while a compact time direction produces a positive Polyakov loop. If more 
than one spatial dimension is short, the bulk volume is reduced and we observe tunneling 
between the two states. The values of the Polyakov loops in the compact direction are not 
correlated, showing that indeed all vacua are equivalent. 

At finite temperature already one spatial dimension allows tunneling of the Polyakov 
line in the compact space dimension. Fig. ^ shows Polyakov lines from such a geometry. 
While those from the non-compact dimensions cluster around zero (P y ; P z not shown) the 
time-like loops have a positive mean value and those in the compact spatial dimensions 
spread along the complex Z{2>) directions. It is not known whether the transition persists 
at high temperature and whether the finite temperature transition at 1/L X = and the 
C breaking transition at zero temperature are connected. However, an exploration of the 
phase diagram of inverse temperature and compact length might be interesting. 
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4. Summary 

For the usual kind of QCD simulations, where the goal is to model QCD in infinite volume, 
the C-breaking transition is just an artifact which needs to be checked. For example, 
Unsal and Yaffe commented that our study of the Nf = 1 condensate might have been 
done in the C-broken phase. We checked our data sets for our studies of the Nf = 1 and 2 
condensates [7, 30] and verified that we were in the C-symmetric phase. 

In the present paper we studied QCD with dynamical quarks with one compact spatial 
dimension. We find a fundamental length scale, the critical length L c , below which QCD 
with four dynamical quark flavors shows spontaneous breaking of charge conjugation. The 
continuum value of L c is approximately 1 fm. It seems to be connected to the finite 
temperature transition in the limit of infinite quark and persists down to the smallest 
quark masses we simulated. There we also observe a crossover in the chiral condensate, 
signaling the restoration at chiral symmetry at the same critical length. We have not 
determined the order of the transition: this would require more extensive finite-size scaling 
studies. 

The spontaneous breaking of charge conjugation is expected from both perturbative 
arguments and a strong coupling expansion. For technical reasons we performed simulations 
with four quark flavors but do not expect the conclusions to change significantly for other 
choices. (We have done some pilot studies with Nf = 2 flavors of improved Wilson fermions, 
and have also seen the crossover to a C-broken phase.) Clearly we are seeing only the tip 
of the iceberg: just as for the finite-temperature transition, the transition is affected by the 
number of colors and flavors, group representation of the fermions, and the relative sizes 
of all small compact dimensions. 
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